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Topological crystalline insulators (TCIs) are insulating materials whose topological property re¬ 
lies on generic crystalline symmetries. Based on first-principles calculations, we study a three- 
dimensional (3D) crystal constructed by stacking two-dimensional TCI layers. Depending on the 
inter-layer interaction, the layered crystal can realize diverse 3D topological phases characterized 
by two mirror Chern numbers (MCNs) defined on inequivalent mirror-invariant planes in 

the Brillouin zone. As an example, we demonstrate that new TCI phases can be realized in layered 
materials such as a PbSe (001) monolayer/h-BN heterostructure and can be tuned by mechanical 
strain. Our results shed light on the role of the MCNs on inequivalent mirror-symmetric planes in 
reciprocal space and open new possibilities for finding new topological materials. 


New topological states of matter, topological crys¬ 
talline insulators (TCIs) PQ, have been identified that 
extend the topological classification beyond the proto¬ 
typical Z2 classification based on time reversal symmetry 
HE]. In TCIs, topological properties of electronic struc¬ 
ture such as the presence of robust metallic surface states 
arise from crystal symmetries instead of time-reversal 
symmetries. There are many proposed TCI phases de¬ 
pending on different crystal symmetries BUS], yet those 
relying on mirror symmetry ns] are of particular inter¬ 
est as they have been experimentally observed in, for 
example, IV-VI semiconductors SnTe, Pbi-^Sn^Te, and 
Pbi-^Sn^Se pTHT 5 ] . More materials are theoretically 
proposed to realize the TCI phases such as rocksalt semi¬ 
conductors da on, pyrochlore iridates PI, graphene 
systems da, heavy fermion compounds [20] [21], and an- 
tiperovskites [ 22 ], including two-dimensional ( 2 D) mate¬ 
rials such as SnTe thin films ^ 3 ^ [ 24 ] and a ( 001 ) mono- 
layer of PbSe [ 25 ] . 

Mirror-symmetric TCIs are mathematically character¬ 
ized by mirror Chern numbers (MCNs). The MCN is a 
topological invariant defined by fii = (/i + —/i_)/2 where 
/i + and /i_ are Chern numbers of Bloch states with the 
opposite eigenvalues of a mirror operator (M z ) calculated 
on the mirror-invariant plane at k z = 0 in the Brillouin 
zone (BZ). In a three-dimensional ( 3 D) crystal, there 
is a second MCN (/i2) defined on the mirror-invariant 
plane at the boundary of the BZ k z = it (in units of 1 /a, 
where a is the length of the primitive lattice vector along 
the z-axis) m- Moreover, considering different mirror 
symmetries, multiple pairs of MCNs (/i 1^2) can be si¬ 
multaneously present in three dimensions. A complete 
characterization of 3 D TCIs requires consideration of all 
the MCNs, which may allow for the possibility of new 
states of matter, where MCNs are locked together or un¬ 
dergo separate transitions. Nonetheless, previous study 
based only on /ii has not explored this situation. 

In this paper, by considering MCNs on all inequivalent 
mirror-symmetric planes in reciprocal space, we study 
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FIG. 1. Schematic drawing of (a) the proposed layered struc¬ 
ture consisting of a stack of 2D TCI layers and (b) the corre¬ 
sponding BZ. Only eight periods of the crystal are shown in 
the crystal structure. The 2D planes at k z = 0 and k z = tt in 
the BZ (gray shaded) are mirror-invariant under the reflection 
, z —> — z, on which the first and second mirror Chern num¬ 
bers (MCNs) are defined, respectively. For a surface normal 
(100), those mirror-invariant planes are projected on Y—F—Y 
and T — Z — T, respectively, along which surface Dirac points 
are expected to occur. 


new topological states of matter realized in a 3 D lay¬ 
ered crystal generated by stacking 2D TCI layers. We 
show that the layered system realizes a new class of 3 D 
TCIs when inter-layer interaction is weak, which we will 
refer to as a layered TCI. The layered TCI is charac¬ 
terized by equal and nonzero first and second MCNs 
Ti = T2 7^ 0 with a number of metallic surface states 
eqaul to |/ii| + |/^2 1 • Increasing the inter-layer interac¬ 
tion, we then show that the layered TCI undergoes topo¬ 
logical phase transitions that change the MCNs (fi 1, M2)- 
Based on first-principles calculations, we predict that a 
heterostructure consisting of alternating layers of PbSe 
monolayer and hexagonal BN (h-BN) sheet realizes the 
layered TCI indexed by (2,2), and that it undergoes 
distinct topological phase transitions in the sequence 
(2,2) —>• (0,2) —)> (0,0) under external uniax¬ 
ial tensile strain. Our findings shed light on new states 
of matter allowed by the presence of multiple MCNs in 
a 3 D crystal. They may also help guide the discovery of 
more topological materials. 

Before presenting the results, we first briefly explain 
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how 2D TCI layers with a non-zero MCN fi^D — n 
(n A 0) can be stacked into a new class of 3D TCIs char¬ 
acterized by /ii = fi^ = n. Consider first a layered system 
consisting of 2D TCIs with /i^d — n 7^ 0 stacked along 
the normal direction to the plane (defined as ^-direction) 
as shown in FIG.[l][a). The layered system then respects 
the mirror symmetry M z that defines the MCN of the 
2D TCI fi^D in the plane of each layer. Now, let us ini¬ 
tially assume that the interaction between the layers is 
negligibly weak, so that every cross section of the 3D BZ 
at constant k z is essentially a copy of the 2D BZ of the 
film. In particular, the mirror invariant planes at k z = 0 
and k z = n [See Fig. [ljb)] should adopt the same MCN 
as the 2D TCI, and thus be indexed by (/i 1 ,^ 2 ) = (n, n). 
For mirror symmetries inequivalent to M z (if any), the 
corresponding MCNs are all trivial (0,0) because the mir¬ 
ror planes allowed by the layered geometry are normal 
to the films, and the crystalline surfaces respecting the 
mirror symmetries are essentially the 2D TCIs without 
metallic (surface) states. This means that the proposed 
TCI are characterized by the coupled MCNs (n, n) for 
M z and (0, 0) for any mirror symmetry inequivalent to 
M z . Turning on the inter-layer interaction in a way that 
respects the mirror symmetry, the MCNs should persist 
within a finite range of the interaction, until the system 
experiences a topological phase transition through a gap 
closure EH, which can lead either to a new topological 
state where the indices are decoupled or to a conventional 
insulating state. 

We demonstrate the topological phases associated with 
MCNs (/ii,/i 2 ) and their transitions from first princi¬ 
ples by applying the above theory to a PbSe/h-BN het¬ 
erostructure. Our calculation is performed with density 
functional theory (DFT) including the Perdew-Burke- 
Ernzerhof [28] generalized gradient approximation as 
implemented in the QUANTUM ESPRESSO package 
[29]. The atomic potentials are modeled by norm- 
conserving, optimized, designed nonlocal pseudopoten¬ 
tials with fully relativistic spin-orbit interaction gener¬ 
ated by the OPIUM package [201 ED- The wave functions 
are expanded in a plane-wave basis with an energy cut¬ 
off of 650 eV. For computational convenience, the energy 
cutoff is reduced to 540 eV when calculating the sur¬ 
face band structure of PbSe(001) monolayers. The van 
der Waals interaction is described based on the semiem- 
pirical dispersion-correction DFT (DFT-D) method [32]. 
The tight-binding model, introduced in Ref. m, is also 
employed to analyze the DFT results on (001) PbSe lay¬ 
ers using parameter sets obtained from our DFT calcula¬ 
tions. The results are consistent with the previous stud¬ 
ies of IV-VI semiconductors, including PbSe TOl l23l [25] . 
The unit cell of the PbSe/h-BN heterostructure is gen¬ 
erated by contracting the in-plane lattice constants of 
the h-BN sheet so that they match the pristine lattice 
constant of the PbSe. We have checked that the artifi¬ 
cial contraction has negligible influence on the electronic 
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FIG. 2. (Color online) Atomic geometry and band structure 
for a slab of PbSe multilayers with a (100) face, (a) The top 
view (upper) and lateral view (left lower) of the slab supercell, 
(b) The band structure for the slab geometry indexed in the 
2D BZ. The MCNs on the k z — 0 and k z — n planes each 
give rise to the 2D Dirac cone on one of the mirror-symmetric 
lines. The gray shaded regions represent the surface-projected 
bulk continuum bands. 

structure near the Fermi energy, as h-BN has a wide band 
gap. 

We first build a layered TCI based on (001) PbSe 
monolayers. Whereas PbSe is a trivial insulator in a 
3D rocksalt geometry, (001) PbSe monolayer is expected 
to be a 2D TCI, indexed by the MCN \/jL 2 d\ = 2 [25] . 
As shown in Fig. (2|a), we consider a system consisting 
of PbSe monolayers stacked along the perpendicular di¬ 
rection to the plane ([001]-direction), so that Pb (and 
Se) atoms form chains along [001], separated by 8.9 A. 
The other crystal parameters are set to those of the bulk 
PbSe. In this way, the inter-layer interaction remains 
weak, and the resulting system is a layered TCI indexed 
by (2, 2) associated with the (001) mirror plane. The sys¬ 
tem respects additional mirror symmetries about {100} 
and {110} mirror planes, on which the MCNs are all triv¬ 
ial as discussed above. 

The calculated MCNs (2, 2) signal the presence of four 
surface states on the mirror-symmetric facets. As de¬ 
picted in Fig. 0 b )> for a surface containing k z , the sur¬ 
face BZ has two inequivalent mirror-symmetric lines, Y 
- T - Y and Z - T - Z which are the projections of the 0 
and 7 r mirror-planes into the surface plane. The absolute 
values of the MCNs |/ii| and |/^ 2 1 dictate the numbers of 
pairs of counter-propagating surface states on the k z = 0 
and k z = tv mirror lines, respectively. It follows that 
there must exist two pairs of surface states along each 
line. To look for the surface states guaranteed by the 
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FIG. 3. (Color online) Topological phase transition in the 
(001) PbSe/h-BN heterostructure. (a) The crystal structure 
and (b) the corresponding BZ. (c) The band structures of the 
PbSe/h-BN heterostructure at h — 3.4 A, h = 4.5 A, h = 2.9 
A, and h = 3.3 A. The MCNs hosted in k z = 0 and k z — ty 
mirror planes at the equilibrium inter-layer distance h = 3.4 
A are adiabatically the same as those at h > 3.4, which is 
a layered TCI with the MCNs (/x 1 , 112 ) = (2,2). The Dirac 
cones at h = 2.9 A, and h — 3.3 A (magnified in the inset) 
respectively signal the topological phase transitions (0,0) —> 
(0,2) and (0,2) (2,2). 


MCNs, we calculate the 2D band structure for the slab 
geometry illustrated in Fig. Ufa). Figure ^b) shows the 
band structure for a slab exposing the (100) surface to 
vacuum along high symmetric lines [See Fig. [ljb)]. As 
expected, in addition to the bulk states in the gray re¬ 
gion, we find surface states that traverse the gap forming 
Dirac points on the mirror-symmetric lines T - Y and T - 
Z. It is clear from the results that the (100) surface has 
four total Dirac points [two shown in Fig. (2|b) and two 
more at the minus of these], dictated by their sum of the 
absolute value \/ii | + |/i 2 1- Note that Y - T and Z -V host 
no metallic surface states. This proves (^ 1 ,^ 2 ) = (0, 0) 
on (010) mirror planes. 

Having demonstrated the layered TCI with fJL\ = fi 2 — 
2 in the layered PbSe system, we now consider a more 
realistic material. Above, we manually fixed the distance 
between the PbSe monolayers to 8.9 A to make the inter¬ 
layer interaction weak, but this is energetically unfavor¬ 
able as the PbSe layers feel a repulsive force as the same 
cations and anions in different layers face each other at 
the interface. To stabilize the system, while keeping the 
inter-layer interaction weak, we insert a h-BN sheet which 
serves as spacer between the neighboring PbSe layers, as 
shown in Fig. [3](a) . A h-BN sheet is a normal insulator 
with a wide band gap of 5 eV, which suggests that the 
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FIG. 4. (Color online) (a) Band gap evolution (lines for 
k z — 0 and k z — tt gap; shading for overall gap) as a function 
of the distance between PbSe and h-BN layers h and (b) corre¬ 
sponding topological phase diagram for the (001) PbSe/h-BN 
heterostructure. ho is the equilibrium position. 


band topology of the heterostructure should be governed 
by bands from the PbSe films. We find that the het¬ 
erostructure has an equilibrium distance (ho) of 3.4 A, 
with a binding energy of 0.08 eV per unit cell of PbSe, 
which indicates that the interaction is in the typical van 
der Waals regime. 

In Fig.[3j we show the band structures of the PbSe/h- 
BN heterostructure along the high-symmetric lines in 
the BZ calculated for various PbSe-h-BN inter-layer dis¬ 
tances (h). First, at equilibrium ho, the system is found 
to be semi-metallic with a small hole pocket at X on the 
k z = 0 plane and an electron pocket near R on the k z = tt 
plane. Then, by increasing the inter-layer distance from 
ho (h > ho), the band gap at X keeps increasing, and the 
system eventually becomes an insulator when h > 3.5 A. 
Increasing h further, we find that the system remains in¬ 
sulating without closing the band gap. Thus, we assign 
the layered TCI phase indexed by (2,2) to the system 
when h > 3.5. 

Conversely, by decreasing the inter-layer distance from 
equilibrium h < ho, which enhances the inter-layer in¬ 
teraction, we find that the system undergoes topological 
phase transitions signaled by the appearance of the Dirac 
points. As presented in Fig. [3jc) , the Dirac points ap¬ 
pear at h = 3.25 A and h = 2.95 A on the k z = 0 
and k z = n mirror planes, respectively. The MCNs, 
calculated using all the valence bands on each mirror- 
symmetric plane, change from (2,2) to (0,2) and from 
(0,2) to (0,0) as h passes through 3.24 A and 2.85 A, re¬ 
spectively. All the topological phase transitions occur in 
a a region where the system is a semimetal because of 
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overlapping bands. Below h = 2.77 A, the valence band 
maximum becomes higher in energy than the conduction 
band minimum on the k z = 0 plane, so the MCN is not 
defined. Therefore, from the strong to weak inter-layer 
interaction regimes, four distinct topological phases ap¬ 
pear, as shown in the phase diagram in Fig. |4^b); a triv¬ 
ial semimetal phase with ( 0 , 0 ), a topological semimetal 
phase with (/i i,/x 2 ) = ( 0 , 2 ), a topological semimetal 
phase with (/i i,// 2 ) = ( 2 , 2 ), and the layered TCI phase 
with (yUi,/i 2 ) = ( 2 , 2 ). Although the heterostructure of 
PbSe/h-BN sheets is expected to be semi-metallic at am¬ 
bient pressure, we expect that these phases should be ac¬ 
cessible under mechanical strain including the proposed 
( 2 , 2 ) layered TCI or by inserting another h-BN sheet be¬ 
tween PbSe layers. We also expect that the phase tran¬ 
sitions demonstrated in this system should be represen¬ 
tative of layered TCIs, and heterostructures of 2 D TCIs 
can be considered as hosts of diverse topological phases 
accessible by engineering the inter-layer interaction. The 
calculated inter-layer distances may vary depending on 
details of the crystal geometry like a stacking registry 
between PbSe-h-BN layers, yet the qualitative features 
should remain intact. 

Finally, we note that layered TCIs are analogous to 
weak topological insulators (TIs) [33] [34]. Weak TIs, 
characterized by zero Z 2 invariant yet non-zero weak 
topological indices (zq zqz/ 3 ), is essentially a stack of 2D 
TI layers along the perpendicular direction that corre¬ 
sponds to G = zqbi + z/ 2 b 2 + z/ 3 b 3 in the BZ 33] . having 
even numbers of robust Dirac cones at the surfaces per¬ 
pendicular to the 2D TI layers [34] 35]. Similarly, 3D 
TCI with the same first and second mirror Chern num¬ 
bers is like layered 2 D TCIs, having |/ii| + |/i 2 | Dirac 
cones on the surfaces normal to the 2 D TCI layers. Also, 
like weak topological indices, (/q, // 2 ) are sensitive to the 
translational symmetry of the crystal. For instance, we 
find that a period-doubling along the z-axis changes the 
MCNs (/ii,/i 2 ): (n, n) —>> (2n, 0), and (rz, — n) —>> ( 0 , 2 n) 
due to the BZ folding, which can be induced by inter¬ 
layer bonding or a registry shift between 2 D TCI layers. 
Indeed, we have found the MCNs to be (4,0) when stack¬ 
ing the PbSe layers with a registry shift atoms between 
Pb and Se atoms in every layer. Notwithstanding the 
sensitivity, it is important to note that the total number 
of surface Dirac cones on the mirror symmetric surfaces, 
dictated by |/ii| + |/^ 21 , is invariant under the period¬ 
doubling, and the second MCN is thus indispensable for 
characterizing the layered TCIs, and indeed all 3D TCIs. 

In conclusion, we propose new topological states of 
matter generated by stacking 2 D TCIs, where simulta¬ 
neous consideration of multiple MCNs is necessary. In 
the non-interacting limit between layers, the layered TCI 
phase emerges where the first and second MCNs (/ii, ^ 2 ) 
are coupled. The layered TCI is a generic class of 3D 
TCIs which can apply to a range of 2 D TCI materials. 
For example, a SnTe thin film, which is expected to be 


a 2 D TCI when cleaved into an odd number of (001) 
layers > 5 [23] . can play the role of the PbSe layer in 
the PbSe/h-BN heterostructure, thus realizing the lay¬ 
ered TCI indexed by ( 2 , 2 ) when the layers are well sep¬ 
arated. Apart from the non-interacting regime, we find 
topological semimetal phases indexed by ( 2 , 2 ) and ( 0 , 2 ) 
and trivial semimetal phase with (0, 0). Despite the pres¬ 
ence of metallic bulk states, the phase transitions should 
be observable via experimental techniques such as angle- 
resolved photoemission spectroscopy. Our findings shed 
light on the possibility of new TCI phases relying on the 
fact that a crystal in three dimensions can have multiple 
MCNs hosted on inequivalent mirror planes in recipro¬ 
cal lattice. These may open the way towards the search 
for new topological materials, based on which quantum 
devices for electronics as well as spintronics can be built. 
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